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AN INTEGRAL TEST OF THE INELASTIC CROSS SECTIONS OF 
Pb AND Mo USING MEASURED NEUTRON SPECTRA 
by D. F. Shook, D. Fieno, C. H. Ford, and G. N. Wrights 

Lewis Research Center 

SUMMARY 

Comparison of measurements and calculations of fast neutron spec- 
tra from a radioactive neutron source inside spheres of Mo or Pb and 
from a cylindrical reactor containing a thick Pb or Mo reflector are 
used as a test of ENDF cross sections. The sphere leakage spectra 
were measured at a sphere -to -spectrometer distance of 2 meters using 
a 54 Ci spherical Arn-Be neutron source. The Pb sphere was 38 cm in 
diameter providing a shell thickness of 4 mfp at 4 MeV. The Mo sphere 
contained powder and was 23 cm in diameter with a shell thickness of 
0. 74 mfp at 4 MeV. Reactor leakage spectrum measurements were 
made at the surface of the NASA ZP-1 reactor which is a cylindrical 
UO 2 F 2 water solution reactor 25 cm in diameter and up to 76 cm high. 
Measurements were made of the bare reactor and with a Pb radial re- 
flector 21 cm thick and with a metallic Mo radial reflector 10 cm thick. 

The leakage spectrum measurements were made in the energy 
range 0. 5 to 11 MeV using a liquid scintillator proton recoil spectrom- 
eter. The spectra were calculated using the S n formulation of the 
neutron transport equation. The flux at the detector was calculated 
using the known source strength in the case of the sphere measurements, 
and per watt of power in the case of the reactor measurements. For Mp, 
ENDF/Bllll microscopic cross sections were used and for Pb, ENDF/ 
B43 cross sections were used. Tests of the adequacy of the scattering 
and angular quadrature order have been made. 

The results for the sphere measurements and calculations show 
that in the case of the thin Mo sphere there is agreement between the 
calculation and measurement. The Pb calculation is much lower than 
the measurement except at the highest neutron energy. 
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In the case of the reactor spectra the use of two-dimensional calcula- 
tions results in spectra that are in good agreement with the measure- 
ments at the surface of the bare reactor indicating that the reactor source 

* 

is reasonably well known. Significant differences in leakage spectrum 
shape are observed for both Mo and Pb reflectors suggesting that there 
are large uncertainties in the inelastic cross sections for Pb and some 
for Mo. 


INTRODUCTION 

Measured and calculated fast neutron spectra in the energy range 
0. 5 to 11 MeV transmitted by lead and molybdenum shields are used as 
a test of existing ENDF/B cross sections for these materials. The 
elements Pb and Mo are tested because of their possible use in compact 
reactors and minimum weight shields for space power applications. Two 
sets of measurements and calculations were made utilizing different 
source spectra. The first set of measurements were made by placing an 
Am-Be (a, n) neutron source inside spheres of the elements tested. 

The sphere leakage spectra were measured at a distance of 2 meters 
using a large liquid scintillator. With this simple geometry it was pos- 
sible to make transport calculations with a large number of energy 
groups and to reduce the uncertainty in the calculation to essentially the 
input cross sections. 

In the second set of measurements the elements tested were used as 
radial reflectors for a small cylindrical homogeneous reactor. The 
spectrum measurements were made close to the outer surface of the re- 
flector to minimize ray effects in the two-dimensional transport calcula- 
tions.' 


Experiment 

The methods used for the sphere measurements are the same as 
reported in reference 1. The sphere leakage spectra were measured at 
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a sphere-to-spectrometer distance of 2 meters with both the sphere and 

spectrometer located 160 cm above the floor of a large room, as shown 

in figure 1- Background due to neutrons scattered into the spectrometer 

from the walls and supporting structure was accounted for using a 

shadow cone. This correction to the measurement was negligible at the 

highest neutron energies measured but increased to 30 percent at the 

lowest energy measured. The measurements were made with a large 

241 

spherical Am -Be neutron source containing 54 Ci of Am (16. 7 g) 

3 

and 66. 8 grams of Be in a net volume of 84. 3 cm . The output of the 

o 

source was (1. 30 ±0. 08)xl0 neutrons per second. The outside diameter 
of the source is 6. 096 cm and includes two stainless steel containment 
shells, one 0. 140-cm thick and the other 0. 152-cm thick. The spherical 
cavity for the source material is 5. 462 cm inside diameter. The Mo 
sphere used in the measurements is 23 cm in diameter and contains Mo 
powder with a density of 3. 73 g/cc providing a shell thickness of 0. 74 mfp 
at 4 MeV. The Pb sphere was 38 cm in diameter providing a shell thick- 
ness of 4 mfp at 4 MeV. The isotopic composition of the Pb in the sphere 
and also the reflector was close to that of natural Pb. The abundances 
are 51. 37 percent Pb 208, 20. 50 percent Pb 207, 26. 83 percent Pb 206, 
and 1. 29 percent Pb 204. 

The neutron spectra were measured using a 5 by 5 cm NE213 liquid 
scintillator proton recoil spectrometer. Pulse shape analysis and two 
parameter data acquisition were used to discriminate against y ray 
pulses. A description of the spectrometer and calibration procedures 
are contained in reference 2. For a scintillator as large as that used 
here, care must be exercised in determining the neutron spectra from 
the measured proton recoil spectra due to multiple scattering and carbon 
interaction in the scintillator. 

In the present work the measured proton recoil spectra were reduced 
using the measurements of Verbinski, et al. (ref. 3) and the unfolding 
code, Ferdor (ref. 4). Since the present measurements are in an es- 
sentially parallel beam geometry, the reference 2 response functions 
are directly applicable. To assure that small differences in scintillator 
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size or spectrometer resolution did not cause a significantly different 
response matrix fit, monoenergetic spectra were measured at 2. 8 and 
14. 7 MeV and detailed comparison satisfactorily made with the refer- 
ence 2 response functions. 

The reactor leakage spectra measurements were made with a 
1. 2 cm diameter by 1. 3 cm long NE218 liquid scintillator positioned with 
the curved surface 3 cm from the reactor surface at the center. The 
neutron spectra were obtained from the measured proton recoil distribu- 
tion using differentiation analysis as outlined in reference 2. 

The reactor, as shown in figure 2 and described in reference 5, was 
cylindrical, 25 cm in diameter and 50 cm high and contained fully en- 
riched uranium as a dilute solution (H/X = 237) of uranyl fluoride in 
water. The reactor was contained in a stainless steel tank having 
0. 159 cm thick walls. Measurements were also made with full radial 
reflectors of lead or molybdenum surrounding the reactor. The metallic 
molybdenum reflector was 10. 1 cm thick and reactor and reflector was 
72 cm high. The lead reflector was 21 cm thick and reactor and reflec- 
tor was 43 cm high. ' v 

The room background was determined to be less than 1 percent of 
the f or egr oun^or • the spectrum measurements at the surface of the re- 
flectors. This was determined from background measurements at 
1 and 2 meters from the reactor utilizing a shadow cone. There is a 
component' of<the measured background due to collisions in the photo- 
multiplier fthbe that would not be independent of source distance, but 
this was: not. evident in the measurements at 1 and 2 meters. Single 
collisions in the photomultiplier tube indicate that this flux is less than 


5 percent of the foreground. 


/ ' TRANSPORT CALCULATIONS 

;V ‘ Sphere Leakage Spectrum 

N 5.V-. 

Multigroup calculations using the S n method were performed to ob- 
tain the sphere leakage spectrum. The calculations were performed in 
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spherical geometry with shells of the test material enclosing the Am-Be 
neutron source. The Am-Be neutron source and its stainless steel con- 
tainers were explicitly included in the calculations. The calculations 
were of the fixed source type using an input spectrum uniformly distrib- 
uted throughout the source region. Fifty-one energy groups were used. 
The first 49 energy groups had a lethargy interval of 0. 1 and covered 
the energy range of 14. 92 to 0. 11 MeV. 

Microscopic cross section data for the test materials (Mo and Pb) 
were obtained from the ENDF/B data files for the group split used. 

The ENDF/B cross sections used for lead were issued as ENDF/B 
material number 43. The data for ENDF/B43 were obtained from 
file 26 of the Aldermaston evaluation. For Mo ENDF/B 1111 cross sec- 
tion were used. The Mo data were originally compiled in 1966 by Pen- 
nington and Gajniak (ref. 6) and issued as ENDF/B material number 
1025. The primary source of information for this evaluation by Pen- 
nington and Gajniak is the evaluation prepared by Schmidt (ref. 7). 

The Mo cross sections have been further revised and issued as 
ENDF/B material number 1111. A number of changes in the ENDF/B 
number 1025 evaluation were made. The most important change for 
the purpose of this study was in the high energy (n, y) cross sections. 

A 3 

These (n, y) cross section data were reduced in value so that computed 
central reactivity worths in fast assemblies correspond to experimental 
reactivity determinations. 

In our calculations the elastic scattering cross sections for Mo in- 
cluded both the Pq and Pj down-scattering transfer cross sections; 
for lead elastic scattering was treated through the Pg order. For Mo 
and Pb the inelastic scattering cross sections, as well as any (n, 2n) 
cross sections, were assumed to have only the Pq component of the 
down- scattering transfer cross sections. 

Calculations for the Mo experiment indicated that a P-^ treatment 
of elastic scattering was adequate for computing the neutron leakage 
spectrum. Calculations in which the elastic scattering is treated 
through the P^ order and through the Pg order differ by less than 
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1. 5 percent for groups 15 through 49 (3. 3 to 0. 11 MeV). The differ- 
ence at group 5 (10 MeV) is about 4 percent. 

Along with the elastic scattering order, the Sn quadrature order 
is important. The calculations reported herein used an S-^g Gauss- 
Legendre quadrature set. This set was found to give an; adequate de- 
scription of the leadage spectrum from the spherical shells of the test 
materials in that a higher order Sgg calculation differs by less than 
1 percent for all groups. 

The calculations reported herein are based on the total number of 

O 

neutrons emitted by the Am -Be source (1. 30 iO. 08x10 neutrons/sec). 

The calculated flux per MeV at 2 meters is folded with the resolu- 
tion function of the spectrometer using a modification of the smoothing 
subroutine in 05S (ref. 8). Thus, the comparison of the calculated 
fluxes with the measured fluxes is on an absolute basis with no nor- 
malization required. 


Reactor Leakage Spectrum 

Multigroup two-dimensional calculations using the Sn method 
were performed to obtain the leakage spectrum for the three critical 
fuel solution reactor configurations considered. The two-dimensional 
calculations were performed in cylindrical (r-z) geometry using DOT 
(ref. 9). Thirty-eight energy groups were used to describe the neutron 
leakage spectrum for the experiments considered. The first 27 groups 
had a lethargy width of 0. 1 and covered the energy range from 14. 92 to 
1. 00 MeV. The next six groups covered the energy range down to 

• 

0. 11 MeV and the remaining groups covered the range of 0. 11 to 0 MeV. 

The calculations for the three critical reactor configurations are 
based on a reactor power of 1 watt. The flux per MeV is calculated at 
the reactor midplane at a radial position of 3. 175 cm (1-| in.) from the 
reactor tank for the unreflected reactor or from the Mo or Pb reflector 
for the reflected reactors. In order to compare the calculations with 
the measured spectra, the resolution function of the spectrometer is 
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folded with these calculated multigroup fluxes as in the case of the 
sphere calculations. 

Microscopic cross section data for the Mo and Pb reflectors were 
obtained from the ENDF/B data files (ENDF/B1111 for Mo, ENDF/B43 
for Pb) while the cross section data for the material components of 
the fuel solution and containment tank were obtained from the GAM-n 
(ref. 10) and GATHER-II (ref. 11) data files. The two-dimensional 
calculations were performed with elastic scattering treated through 
Pg and the angular quadrature through Sg. 

RESULTS AND DISCUSSION 
Bare Source 

The measured and calculated (a, n) source spectrum is shown in 
figure 3 and has been discussed in reference 1. A measured spectra 
for a small (a, n) source was used as input for transmission through 
the larger source. The comparison is generally good but the calcula- 
tion has a few percent less total flux above 7 MeV and smoothes out the 
measured shape changes in this energy region. The {a, n) spectra is 
not expected to have abrupt changes in shape in this region but it does 
have more structure than can be represented by the particular multi- 
group structure used. At energies less than 7 MeV the agreement be- 
tween the calculation and measurement is quite satisfactory. 


Mo Sphere 

The measured and calculated leakage spectrum from a thin Mo 
sphere containing the (a, n) source is shown in figure 3. The agree- 
ment between experiment and calculation is in general satisfactory with 
the calculation slightly higher than the measurement around 1. 0 MeV 
neutron energy. Shape differences between measured and calculated 
leakage spectra for a thick molybdenum reflector on the solution reac- 
tor were observed below 5 MeV and are shown in figure 7. 
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Pb Sphere 

The measured and calculated leakage spectrum from the natural 
lead sphere is shown in figure 3. A comparison between the two spectra 
indicates that the calculation is generally lower than the measurement. 
The agreement is particularly poor at 2. 0 MeV where the calculated flux 
is only about 50 percent of the measured flux. Since the flux at this low 
energy depends not only on the inelastic removal but on the inelastic 
transfer cross sections from higher energies, it is difficult to determine 
the origin of such a large discrepancy from this measurement. A spec- 
trum measurement and calculation made at ORNL (ref. 12) also shows a 
significant undercalculation of the flux around 2 MeV. Since the ORNL 
results depend only on the total cross section,, they indicate that a part 
of the undercalculation around 2 MeV for this experiment is because the 
cross sections are too high in this energy range. 

Bare Reactor Leakage Spectrum 

The interpretation of the Pb and Mo reflector leakage spectrum is 

more complex than the sphere spectra because of the two-dimensional 

geometry and the presence of the reactor core materials of hydrogen, 

235 

oxygen and small amounts of U and fluorine. To aid in the inter- 
pretation a measurement and calculation of the bare core leakage spec- 
tra was made. Tne measured and calculated spectra are shown in fig- 
ure 4. The agreement between the experiment and the DOT calculation 
is seen to be generally satisfactory with the calculation slightly lower 
than the measurement. 

Mo Reflected Reactor Leakage Spectrum 

The measured and calculated leakage spectrum at 3 cm from the 
surface of the molybdenum reflected reactor is shown in figure 5. A 
comparison of the results in the figure show good agreement between 
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the transport calculations and measurement at energies greater than 
6 MeV but significant shape difference at energies less than 6 MeV. 
These shape differences were not apparent in the sphere experiment so 
they could be due in part to processes occurring above a neutron energy 
of 10 MeV where there are no neutrons in the (a, n) source spectra. 
An additional factor is that the Mo reflector is much thicker than the Mo 
shell used in the sphere measurements. 

Pb Reflected Reactor Leakage Spectrum 

The measured and calculated leakage spectrum at 3 cm from the 
surface of the natural Pb reflected reactor is shown in figure 6. A 
comparison of the results in the figure show a discrepancy exists be- 
tween the transport calculation and the measurement that is similar to 
the sphere measurement (see fig. 3). The large discrepancy in the 
energy range 1 to 2 MeV evident in the sphere measurement is, also 
present in this measurement. 
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Figure 3. - Leakage spectrum for Am Be source enclosed 
by a Pb and a Mo shell. 



Figure 4. - Leakage spectrum at surface of bare cylindrical 
solution reactor. 




NEUTRON ENERGY, MeV 

Figure 5. - Leakage spectrum at surface of molybdenum 
reflected reactor. 



figure 6. - Leakage spectrum at surface of lead reflected 
reactor. 
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